# Required fxn: 'cluster1 n'

# elements: mi, income (response): yi, rent (y/n): ai
M=2500; N=415
s2r=25189^2
s2t=21784^2
s2p=.726^2


# find n for mean with known M
Mbar=M/N
c1.n('income',s2r,Mbar,N,'total',1000000)
#  or estimate n with estimated Mbar
Mbar=6.04
c1.n('income',s2r,Mbar,N,'total',1000000)

# estimate n with unknown M
c1.n('income',s2t,6.04,N,'total',1000000)

# estimate proportion
c1.n('income',s2p,Mbar,N,'proportion',0.04)